Effects of macroscopic polarization in III-V nitride multi-quantum-wells 
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Huge built-in electric fields have been predicted to exist in wurtzite III-V nitrides thin films and 
multilayers. Such fields originate from heterointerface discontinuities of the macroscopic bulk polar- 
ization of the nitrides. Here we discuss the background theory, the role of spontaneous polarization 
in this context, and the practical implications of built-in polarization fields in nitride nanostructures. 
To support our arguments, we present detailed self-consistent tight-binding simulations of typical 
nitride QW structures in which polarization effects are dominant. 
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I. INTRODUCTION 

Spontaneous polarization has long been known to take 
place in ferroelectrics. On the other hand, its existence 
in semiconductors with sufficiently low crystal symme- 
try (wurtzite, at the very least) has been generally re- 
garded as of purely theoretical interest. Recently, a se- 
ries of first principles calculations!!! - !] has reopened this 
issue for the technologically relevant III-V nitride semi- 
conductors, whosa natural crystal structure is, in fact, 
wurtzite. Firstly,EI it was shown that the nitrides have 
a very large spontaneous polarization, as rWjU as large 
piezoelectric coupling constants. Secondlyp'El it was di- 
rectly demonstrated how polarization actually manifests 
itself as electrostatic fields in nitride multilayers, due to 
the polarization charges resulting from polarization dis- 
continuities at heterointerfaces. This charge-polarization 
relation, counterchecked in actual ab initio calculatians,tl 
has been exploited to calculate dielectric constants.13 

While piezoelectricity-related properties are largely 
standard, spontaneous polarization is to some extent 
new in semiconductor physics, to the point that, so far, 
the practical importance of spontaneous polarization in 
III-V nitrides nanostructures (multi quantum wells, or 
MQWs, are the focus of this paper) has been largely 
overlooked. It is tantalizingly clear to us, however, that 
these concepts may lead to a direct and unambiguous 
measurement of the spontaneous polarization in semicon- 
ductors, to the recognition of its importance in nitride- 
based nanostructures, and, hopefully, to its exploitation 
in device applications. Our work has already spawned 
interpretations] and purposely planned experimentsu in 
this direction; in the hope of accelerating the process, 
in this paper we show how to account for the effects of 
spontaneous polarization in MQWs, and discuss some 



prototypical cases and their possible experimental real- 
ization. To support our arguments, we present simula- 
tions of typical AlGaN/GaN MQWs where spontaneous 
and piezoelectric polarizations are about equal. 

Among the consequences of macroscopic polarization 
which we will demonstrate in this paper let us mention 
the following: (a) the field caused by the fixed polar- 
ization charge, superimposed on the compositional con- 
finement potential of the MQW, red-shifts dramatically 
the transition energies and strongly suppresses interband 
transitions as the well thickness increases; (b) the ef- 
fects of thermal carrier screening are negligible in typical 
MQWs, although not in massive samples; (c) a quasi- 
flat-band MQW profile can be approximately recovered 
(i.e. polarization fields can be screened) only in the pres- 
ence of very high free-carrier densities, appreciably larger 
than those typical of semiconductor laser structures; (d) 
even in the latter case, transition probabilities remain 
considerably smaller than the ideal flat band values, and 
this reduces quantum efficiency; (e) once an appropri- 
ate screening density (i.e. the pumping power or injec- 
tion current) has been chosen to ensure that the recom- 
bination rate is sufficient, a residual polarization fields 
is typically still present: this provides a means to in- 
tentionally red-shift transition energies by changing well 
thicknesses, without changing composition; (f) the very 
existence of distinct and separately controllable sponta- 
neous and piezoelectric polarization components allows 
to choose a composition such that they cancel each other 
out, leading to flat-band conditions. Analogously, for a 
proper choice of superlattice composition, piezoelectric 
polarization can be made to vanish and hence a measure 
of spontaneous polarization can be accessed, through e.g. 
the changes in optical spectra. It is clear that a fuller un- 
derstanding of these points will ultimately lead both to 
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improvements in design and operation of real nitride de- 
vices, and to the direct measurement of polarization, and 
a better knowledge thereof, in nitride semiconductors. 

Before moving on, let us mention other recent contri- 
butions in this area. Buongiorno Nardelli, Rapcewicz, 
and Bernholc, using non-self-consistent effective- mass 
based perturbation theory in the small field limit, have 
predictedB red shifts and transition probability suppres- 
sion in InGaN quantum well; Delia Sala et ai, using 
self-consistent tight binding calculations, have appliedO 
some of the ideas reported in this paper to InGaN/GaN 
quantum well lasers, explaining several puzzling exper- 
imental features, among which the high thresholds ob- 
served for GaN-based lasers, and several other aspects 
related to .self-consistent screening effects; Montecarlo 
simulationsa by Oberhuber, Vogl, and Zandler, employ- 
ing the polarization calculated in Ref. [l[ have revealed 
a polarization-enhanced carrier density in the conduc- 
tion channel of _AlGaN/GaN HEMTs. Takeuchi et al. 
have interpretedu their own measurement of a quantum- 
confined Stark effect in InGaN/GaN superlattices as 
piezoelectricity-induced; this is, as it turns out, essen- 
tially correct in their specific case involving InGaN, but 
not in general. Analogous results have been reported 
in Refs. and [ll], the latter including fairly detailed 
simulations, but only accounting for piezoelectricity. Fi- 
nally, a detailed theoretical exposition of effective-mass 
theory adapted to deal with piezoelectric fields is given 
in Ref. including useful notation and basic formu- 
las, and some applications. Experimental work will be 
mentioned later; here let us just quote the very recent 
evidence that polarization-related effects on optical prop- 
erties in selected AlGaN/GaN systems cannot be prop- 
erly interpreted if spontaneous polarization is neglected; 
circumstantial evidence was obtainedu by Leroux et ai, 
while a more carefully planned investigation, reaching 
firmer conclusions, has been carried outEl by Cingolani et 
al. 



II. PIEZOELECTRIC FIELDS 

Piezoelectricity is a well known concept in semiconduc- 
tor physics. Binary compounds of strategic technologi- 
cal importance as III-V arsenides and phosphides can be 
forced to exhibit piezoelectric polarization fields by im- 
posing upon them a strain field. 

Among others, applications of piezoelectric effects in 
semiconductor nanotechnology exist in the area of multi 
quantum- well (MQW) devices. A thin semicondutor 
layer (active layer) is embedded in a semiconductor ma- 
trix (cladding layers) having a different lattice constant. 
If pseudomophic growth occurs, the active layer will be 
strained and therefore subjected to a piezoelectric polar- 
ization field, which can be computed as 

P (pz) = V • e (1) 



if the strain field e and the piezoelectric constants tensor 
e are known. 

In a finite system, the existence of a polarization field 
implies the presence of electric fields. For the piezoelec- 
tric case, the magnitude of the latter depends on strain, 
piezoelectric constants , and (crucially) on device geom- 
etry. The structure of a typical III-V nitride-based su- 
perlattice or MQW is -C-A-C-A-C-A-C- (A=active, 
C=cladding), where both the cladding layer and the ac- 
tive layer are in general strained to comply with the sub- 
strate in-plane lattice parameter. In such a structure, the 
electric fields in the A and C layers are 

Ei pz) = 4^ c (P<? z) - Pf } ) /{i c e A + £ A s c ) 

E £ z) = 4^ A (Pi pz) - P[? z) ) /{l c e A + £ A s c ) (2) 

where £a.c ar e the dielectric constants and £ A .c the 
thicknesses of layers A and C. Thus, in general, an elec- 
tric field will be present whenewar Pa =/= Pc- The above 
expressions are easily obtainedEj by the conditions that 
the electric displacement be conserved along the growth 
axis, and by the boundary conditions that the potential 
energy on thc_far right and left of the MQW structure 
are the same£J 

There are essentially three special cases of MQW struc- 
tures worth mentioning: 

i) active (cladding) layer lattice matched to the sub- 
strate: P A = (P c = 0); 

ii) £ A — £q, whence Ea = — Eq; 

iii) £ A <C £c, implying Ec — 0, and hence 

Ei pz) = 47rPi pz) /£A. (3) 

In the last case we implicitly assumed the cladding layer 
to be unstrained - that is, its lattice constant to be re- 
laxed to its equilibrium value because its thickness ex- 
ceeds the critical value for pseudomorphic growth over 
the substrate. p( pz ) may take any direction in general, 
but in norrua.1 practice its direction is parallel to the 
growth axis.li3 

To obtain piezoelectric polarization effects in 
zincblende semiconductor systems, lattice-mismatched 
epitaxial layers are purposely grown along a polar axis, 
e.g. (HI); the in-plane strain propagates elastically 
onto the growth direction, thereby generating p(? z ). 
In wurtzite nitrides, the preferred growth direction is 
the polar (0001) [or (0001)] axis, so that any non- 
accomodated in-plane mismatch automatically generates 
a piezoelectric polarization along the growth axis (the 
sign depends on whether the epitaxial strain is compres- 
sive or tensile). We will be always be assuming this 
situation in the following. 

Usually, alloys are employed in the fabrication of 
MQWs. In that case, one may estimate the piezoelec- 
tric polarization in the spirit of Vegard's law as, for a 
general strain imposed upon e.g. an Al^Gai-^N alloy, 
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P( pz > = [x e*AiN + (l-x) e GaN ] e{x) , (4) 

This expression contains terms linear as well as quadratic 
in x, and similar relations hold for quaternary alloys. 
This piezoelectric term is only present in pseudomorphic 
strained growth, and will typically tend to zero beyond 
the critical thickness at which stain relaxation sets in. 
Uncomfortable though it may be,t3 the Vegard hypotesis 
is at this point in time the only way to account for piezo- 
electric (and spontaneous, see below) fields in alloys. As 
will be shown below, indeed, the qualitative picture does 
not depend so much on the detailed value of the polar- 
ization field as on their order of magnitude. 



III. SPONTANEOUS FIELDS IN MQWS 

New possibilities are opened by the use of III-V nitrides 
(AlN,GaN,InN), that naturally cristallize in the wurtzitc 
structure. These materials are characterized by polariza- 
tion properties that differ dramatically from those of the 
standard III-V compounds, considered so far. From sim- 
ple symmetry arguments,li3 it can be shown that wurtzite 
semiconductors are characterized by a non-zero polariza- 
tion in their equilibrium (unstrained) geometry, named 
spontaneous polarization (or, occasionally, pyrpelectric, 
with reference to its change with temperature) J13 While 
the spontaneous polarization of ferroelectrics can be mea- 
sured via an hysteresis cycle, in a wurtzite this cannot 
be done, since no hysteresis can take place in that struc- 
ture. Indeed, spontaneous polarization has never been 
measured directly in bulk wurtzites so far. III-V ni- 
tride MQWs offer the opportunity to reveal its existence 
and to actually measure it. In turn, spontaneous po- 
larization can provide new degrees of freedom, in the 
form of permanent strain-independent built-in electro- 
static fields, to tailor transport and optical character- 
istics of nitride nanostructures. Its presence can e.g. be 
exploited to cancel out the piezoelectric fields produced 
in typical strained nitride structures, as discussed below. 

Thanks to recent advancesEj in the modern theory of 
polarization (a unified approach based on the Berry's 
phase concept), it has become possible to compute eas- 
ily and accurately from first principles the values of the 
spontaneous polarization, besides piezoelectric and di- 
electric constants, in III-V nitrides .Em The results of 
the calculations show that III-V nitrides have important 
polarization-related properties that set them apart from 
standard zincblende III-V semiconductors: 

i) huge piezoelectric constants (much larger than, and 
opposite in sign to those of all other III-V's); 

ii) existence of a spontaneous polarization of the same 
order of magnitude as in ferroelectrics. 

The latter is, we think, a most relevant property. Spon- 
taneous polarization implies that even in heterostructure 



systems where active and cladding layers are both lattice- 
matched to the substrate (so that no strain occurs, hence 
no piezoelectricity), an electric field will nevertheless ex- 
ist due to spontaneous polarization. In addition, unlike 
piezoelectric polarization, spontaneous polarization has a 
fixed direction in the crystal: in wurtzites it is the (0001) 
axis, which is (as mentioned previously) the growth di- 
rection of choice for nitrides epitaxy. Therefore the field 
resulting from spontaneous polarization will point along 
the growth direction, and this (a) maximizes spontaneous 
polarization effects in these systems, and (b) renders the 
problem effectively one-dimensional. In the simplest case 
of a fully unstrained (substrate lattice-matched) MQW, 
the electric fields inside the layers are given, in analogy 
to Eq. | by 

Ei sp) = 4^c(P^ sp) ~ Pi sp) )/(^A + tKSc) 

E^ sp) = 4^ A (pf > - Pg p) )/(^A + Ias c ) (5) 

where the superscript (sp) stands for spontaneous; typi- 
cal spontaneous polarization valuesEJ indicate that these 
fields are very large (up to several MV/cm). 

In actual applications (for instance, to produce un- 
strained MQWs) alloys will have to be employed. The 
values of the spontaneous polarization are accurately 
known only for binary compounds.!!! In the absence of bet- 
ter estimates, we assume as before that the spontaneous 
polarization in alloys can be estimated using a Vegard- 
like rule as (for, e.g., Al x l&yG&i- x -yN) 

Pta0 (x, y) = x P^ + V +(!"*-») Pg& ■ 

In Figure |l] we report the resulting spontaneous polariza- 
tion vs. lattice constant for the III-V nitrides, with data 
from Ref. 
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FIG. 1. Spontaneous polarization in A\ x \riyG&\- x -y^ al- 
loys according to a Vegard- like rule. 

Figure |l| shows that for a given (substrate) lattice con- 
stant, a wide interval of spontaneous polarizations (hence 
of spontaneous fields, according to Eq. ||) is accessible 
varying the alloy composition. In particular, consider 
a GaN/ALJn^Gai-z-yN MQW, where the composition 
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is chosen so that the alloy be lattice matched to GaN, 
which we assume to be also the substrate (or buffer) ma- 
terial (dashed-dotted line in Figure Q). Then, piezoelec- 
tric polarization vanishes, but spontaneous polarization 
remains, and takes on values up to ^0.05 C/m 2 . For 
a GaN quantum well with thick AlGaN cladding layers, 
this means a theoretical electrostatic field of up to about 
5 MV/cm. 



Note again, at his point, that the fields (see Eqs. [2] and 
||) are not related to just the polarization of the material 
composing the specific layer, but a combination of polar- 
ization ^differences, dielectric screening, and geometrical 
factorsEy We now consider the field values in the active 
layer: the total field Ea is shown in Figure || vs. Al and 
In molar fraction; the same is done for the piezoelectric 
component in Figure ||. In both cases the appropriate 
Vegard-like rules have been used. 



IV. FIELDS IN THE GENERAL CASE 

In general, of course, MQWs will be strained. Then, 
for an arbitrary strain state, the electric fields in the A 
(or C) layers of the MQW are the sum of the piezoelectric 
and spontaneous contributions: 
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where E^ pz ^ is the old-fashioned piezoelectric field from 
Eq. ||, and E^ sp ) is given by Eq. ||. It is important to 
stress that E^ sp ^ depends only on material composition 
and not on the strain state. Also, it is a key point to 
notice that although both polarization contributions lay 
along the same direction [the (0001) axis], p( pz ) may have 
(due to its strain dependence) the same or the opposite 
sign with respect to the fixed p( sp ) depending on the 
epitaxial relations. 
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FIG. 2. Total built-in electrostatic field in the active layer 
of a AblnyGai-^-aN/GaN MQW system (see text) vs. Al 
and In molar fraction. Fields are in units of MV/cm, positive 
fields point in the (0001) direction (Ga-face). 

It is difficult to give a simple picture of the electric 
field pattern in a general MQW system because of the 
many degrees of freedom involved. Here we consider an 
Al^Ga^Ini-x-yN/GaN MQW pseudomorphically grown 
over a GaN substrate, having active and cladding layers 
of the same thickness. In such a case 
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FIG. 3. Piezoelectric component of the electrostatic fields 
in the same MQW system as in FigureH (see text). 



Comparison of these Figures cleary bears out the 
importance of spontaneous polarization in determining 
the electric field. Several aspects are worth pointing 
out. First, large electric fields (~ 0.5-1 MV/cm) can 
be obtained already for modest Al and In concentra- 
tions. Second, it is easy to access compositions such that 
Al x In a Gai_2;_ y N is lattice matched to GaN: thereby, no 
piezoelectric fields exist, but large, purely spontaneous 
fields still do; specifically, this situation is realized for 
compositions laying on the zero-piezoelectric-field line in 
Figure |[ On this locus of compositions, spontaneous po- 
larization is the only source of field and it can therefore be 
measured via the changes it induces in the MQW spectra. 
Third, it is possible to choose the material composition 
in such a way that the active layers of a MQWs are free 
of electric fields. To achieve this situation the MQW 
must be strained so that the piezoelectric and sponta- 
neous polarizations cancel each other out; clearly, this 
is realized for compositions laying on the zero-field line 
in Figure ||. Of course the possibility of having a null 
field where desired is of capital importance in those de- 
vices where electric fields in the active layer can not be 
tolerated (other field screening mechanisms are discussed 
below). 

A noticeable feature of Figure || is that the piezoelec- 
tric component increases much faster with In content that 
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with Al content, despite the larger piezoelectric constants 
of the latter. The reason is, of course, that strain builds 
up much more rapidly with In concentration. Along with 
the small difference in spontaneous polarization between 
InN and GaN, this is the reason why it is possible to in- 
terpret with reasonable accuracy polarization effects in 
InGaN/GaN structures on the basis of purely piezoelec- 
tric effects, as done in Refs. H-JLl]. On the other hand, 
it can be seen that the spontaneous component increases 
much more rapidly with Al content than with In con- 
tent, due to the widely different polarizations of A1N and 
GaN. For AlGaN, piezoelectricity-based interpretations 
are bound to fail. 



V. EFFECTS OF POLARIZATION FIELDS 

We now come to the implications of polarization fields 
for devices based on III-V nitrides. In this Section we 
present a set of accurate self-consistent tight-binding cal- 
culations for an isolated AlGaN/GaN QW representing 
a system in which the the spontaneous-polarization con- 
tribution to the total built-in electrostatic field is as 
large as the piezoelectric term. To simulate realistically 
these nanostructures, self-consistency is needed to de- 
scribe field screening by free carriers; the latter cannot 
physically cancel out the polarization charge, which is 
fixed and invariable, but may screen it out in part. In our 
calculations we therefore solve self-consistently the Pois- 
son equation and the Schrodinger equation for a state- 
of-the-art empiricaLtight binding Hamiltonian for realis- 
tic nanostructures t3 In the following, two cases are con- 
sidered: (a) non-equibrium carrier distribution (Subsec. 
A and B) related to photoexcitation or injection, where 
electron and hole quasi-Fermi levels are calculated for 
a given areal charge density (u2d) in the quantum well 
(the sheet density, related to the injection current or opti- 
cal pumping power); (b) thermal equilibrium distribution 
(Subsec. C and D) where the Fermi level is calculated as 
a function of doping density by imposing charge neutral- 
ity conditions^ We solve Poisson's equation, 
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where the (position-dependent) quantities D 7 e, and V, 
are respectively the displacement field, dielectric con- 
stant, and potential. Pr is the (position-dependent) total 
transverse polarization. The effects of composition, po- 
larization, and free carrier screening are thus included 
in full. [Consistently with the aim of describing a single 
QW, we choose boundary conditions of zero field at the 
ends of the simulation region. This corresponds to the 
Iq — > oo limit in Eqs. || and ||.] 

The potential thus obtained is inserted in the 
Schrodinger equation, which is solved diagcoializing an 
empirical tight-binding sp 3 d 5 s* Hamiltonian. The pro- 
cedure is iterated to self-consistency. Further appli- 



cations and— details on the technique can be found 
elsewhere. QO 

Here we concentrate in particular on the polarization- 
induced quantum-confined Stark effect (QCSE) in zero 
external field, its control and quenching, and its evolu- 
tion with layer thickness. We first deal with the low 
free-carrier densities regime: thereby the QCSE mani- 
fests itself as a strong red shift of the interband transi- 
tion energy, with a concurrent suppression of the transi- 
tion probability, both of these features getting stronger 
as the well thickness increases. This is the regime that 
applies to low-power operation or photoluminescence ex- 
periments. 

Next we discuss how the QCSE can be modified, and 
eventually (almost) quenched, by providing the QW with 
a sufficiently high free-carrier density. In this regime, 
as the free carrier density increases, the transition en- 
ergy is progressively blue-shifted back towards its flat 
band value, and the transition probability suppression 
is largely removed. The needed free-carrier density de- 
pends on the polarization field, and not surprisingly it is 
found to be very substantial. Typical values of the sheet 
density range in the 10 13 cm -2 , as opposed to typical 
values of 10 12 cm' 2 needed to obtain lasing in GaAs-like 
materials. 



A. QCSE at low power 



The prototypical system we consider is an isolated GaN 
quantum well cladded between Al^Ga^^N barriers. In 
Figure [I], we display the total field Ea in the (isolated) 
active well, and its piezoelectric component as a function 
of the Al molar fraction x. The spontaneous component 
is the difference of the two, and. therefore approximately 
equal to the piezoelectric oneO 

The value we pick for our simulations is x—0.2, a rea- 
sonable compromise between the conflicting needs for 
not-too-large fields, sufficient confinement, and techno- 
logically achievable composition. In this case the valence 
offset is AE V = 0.064 eV. The total field in the QW 
of -2.26 MV/cm, and the spontaneous and piezoelectric 
components are -1.14 MV/cm and -1.12 MV/cm respec- 
tively. The minus signs indicates that the field points 
in the (0001) direction. The bare polarization charge at 
the interface is proportional to the change in polariza- 
tion across the interfaces, and it amounts to ~1.28xl0 13 
cm~ 2 . The field value mentioned above results from this 
charge as screened by the dielectric response of the QW 
(the field change at the interface is thus related to a 
smaller, or screened, effective interface charged). 
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have the character of a resonance for small x {i.e. small 
confinement). 
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FIG. 4. Total field and its piezoelectric component in the 
GaN/ALcGai-zN QW discussed in the text. 

We performed a series of calculations for different well 
widths, where the electron and hole confined states have 
been populated {i.e. pairs have been created) with a den- 
sity of ~ 10 11 cm~ 2 to simulate a low-power optical exci- 
tation. We find that this density has only a very marginal 
effect: indeed, the potential is perfectly linear, i.e. the 
electrostatic field remains uniform, over the whole QW. 
The square-to-triangular change in the potential shape 
causes a small blue shift of both the electron and hole 
confined states (referred to the flat well bottom), but the 
linear potential given by the field causes a much larger 
relative red shift for any reasonable thickness. Also, since 
the thermal carrier density fluctuations are negligible at 
microscopic thicknesses and room temperature (see be- 
low, and Ref. |24|) , one expects the QW band edge profile 
to remain linear as a function of thickness, at least for 
the low excitation powers typical of photoluminescence 
spectroscopy. 

In Figure || we show the TB result for the lowest in- 
terband transition energy and the corresponding transi- 
tion probability {i.e. the squared overlap of the highest 
hole leveL-and the lowest electron level envelope wave- 
functionscil) as a function of QW thickness. Both the 
Stark red shift and the strong suppression of the tran- 
sition probability are evident. This was to be expected 
from the potential shape and the reduced overlap of hole 
and electron states, displayed in the inset of Figure ||. 

It is worth noting that the localization of the hole en- 
velope function in the well region is rather weak, because 
the large effective field blue-shifts the hole bound state 
energy close to the valence barrier edge. This will gener- 
ally be the case for low-cc AlGaN wells, due to the small 
valence confinement energyB Indeed, even the conduction 
confinement is small on the scale of the fields- induced po- 
tential drop, and the electron bound state also tends to 
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FIG. 5. Transition energy red shift and suppression of tran- 
sition probability vs well thickness. In the inset self-consistent 
band edge (solid), electron and hole envelope functions of the 
TB wavefunctions (dashed) for a 50 A thick QW. 



We conclude that in the absence of excitation and at 
normal operation temperatures, or at low optical excita- 
tion power, macroscopic polarization fields cause QW's 
to be highly inefficient in emitting light, and the emission 
energy to be considerably different than the gap of the 
material plus the confinement energies. 

Comparison with experiment is tricky since most at- 
tempts to measure these effects are polluted by inappro- 
priate (at least for the purpose of revealing polarization 
effects) choices of the experimental geometry. For in- 
stance, measurements have been donea in a series of 
quantum wells_pf different thicknesses ranging between 
10 and^Axa In any case, the general experimental 
featuresullj'BEj are in full agreement with the notion that 
the transitions are red-shifted essentially linearly with in- 
creasing well thickness, and that screening at low free car- 
rier densities is irrelevant in this class of systems. This 
is not quite tru e an y more for thick layers, as will be 
discussed in Sec.VC. 
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B. QCSE quenching at high excitation power 

If carriers are generated optically, one can envisage 
that a sufficiently high excitation power could possibly 
produce the carrier density needed to screen the polar- 
ization held. We now calculate the properties of the QW 
as a function of the free-carrier areal density, to check if 
the red shift and the transition probability suppression 
can be removed in a physically accessible range of such 
density. 

We repeat the self-consistent procedure increasing pro- 
gressively the free charge density in the QW from 10 12 
up to 2 x 10 13 cm -2 . We see in Figure [| that, albeit at 
the cost of a large increase of the QW free-carrier density, 
the held does get progressively screened. 

4.5 I ' 1 , 1 | 
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FIG. 6. Self-consistent band edge (a) and electric field (b) 
for various sheet densities (~ excitation level) in a 50 A thick 
QW. 

As can be seen in Figure [j], at fixed thickness the red 
shift decreases as the carrier density increases, and it 
tends to become thickness-independent at the highest 
densities. The transition probability is also increased by 
several orders of magnitude. However, the field is not 
screened abruptly but dies off gradually, with an effec- 
tive screening length of about 20 A for the largest density 
used here (of course, this is a token of the larger spatial 
extension of the screening charge as compared to the po- 
larization charged): herefore, holes and electron remain 
spatially separated to a large extent even at high carrier 
densities, and the trasition probability never quite goes 
back to unity. This is likely to be one of the reasons for 
the relatively low quantum efficiency observed in typical 
nitride MQW devices. For the same reasons, the transi- 



tion energy never quite goes back to the flat-band value 
(gap plus confinement energy). Note in passing that be- 
cause of strain, in these calculations i?^ aN = 3.71 eV, 
almost 10 % larger than the equilibrium value. 
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FIG. 7. Removal of red shift and recovery of interband 
transition probability upon high excitation. 

The screening density of order 2 xlO 13 cm~ 2 needed 
to partially screen out the field corresponds to an 
estimatedEO optical pumping power of about 10 to 20 
kW/cm 2 per well. This figure -agrees mcely with the un- 
usually high pumping powersQ'E£rE3o needed to obtain 
the laser effect in nitride structures. The explanation is 
simply that much of the free charge being generated ac- 
tually goes into screening the polarization field. On the 
other hand, our result prove that the optically activated 
lasing conditions can indeed be realized in practice, al- 
though with high pumping powers, so that thepe seems to 
be no need to invoke quantum dot formations or other 
exotic effects to explain lasing in nitride structures. On 
the other hand, the same phenomenon explains the high 
current threshold observed for electrically driven GaN 
based lasers. EjTJ 

QCSE quenching phenomena similar to those just de- 
scribed have been observedtj by Takeuchi et al. in In- 
GaN/Gan MQWs, with estimated fields in the 1 MV/cm 
range. The red shift and optical inefficiency were in 
fact removed , although only in part and in a tran- 
sient fashion, by sufficiently high excitation powers. The 
order of magnitude of the values reported in Ref. [| is 
~200 kW/cm 2 for 5 to 10 MQW periods, i.e. 20 to 40 
kW/cm 2 per well, in qualitative agreement with our es- 
timate above. 

One important remark at this point is that, depending 
on the excitation power, the MQW will adsorb radia- 
tion at many different transition energies ranging from 
that of the built-in-field-biased well (low power limit) to 
the quasi-flat-band well (high power limit) - that is, the 
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MQW acts a multistable switch. It is indeed fortunate 
that the typical fields in these structures are such that 
one can physically access the various possible regimes. 

Another noticeable effect is that at a properly chosen 
value of the sheet density (i.e. of the excitation power) 
one can obtain at the same time a reasonable transition 
probability and a red-shifted energy by just increasing 
the well thickness. This is very useful since the transi- 
tion wavelength can be shifted to a different color without 
changing alloy composition, but only the well thickness. 
For instance (see Figure changing the well thickness 
from 20 to 30 A at a sheet density of 4xl0 12 cm~ 2 , one 
obtains an energy red shift of 0.1 eV at the cost of a loss 
of a factor 10 in recombination rate, which may still be 
acceptable depending on the application. Red-shifting 
the transition energy in this fashion may avoid the need 
to add e.g. some In in the QW composition. Of course, 
Wite-shifting by thickness reduction will increase the tran- 
sition probability. 

C. Self-screening of fields in massive samples 

Free charge produced by high excitation screens po- 
larization fields fairly efficiently over the quite short dis- 
tances typical in nanostructures, basically because the 
spatial extension of the screening charge is comparable 
to the size of the system. How about extended samples, 
especially if not subjected to illumination, i.e. having 
only intrinsic free carriers ? It is indeed the case that no 
macroscopic fields exist in "infinitely large" samples even 
in the absence of high densities of (say) photogenerated 
carriers. The simple reason is that the intrinsic carrier 
fluctuations in an undoped semiconductor rise exponen- 
tially as a function of deviations of the chemical poten- 
tial from the mid-gap value.L3 In polarized nitrides, such 
deviations occur due to the built-in fields. As the sam- 
ple thickness increases, the potential drop grows linearly. 
When the drop is smaller than the gap, the field is uni- 
form: |E| = 47rP/eo. When the drop approaches the gap 
value, i.e. for thicknesses approaching d c = £g ap /|E|, 
the Fermi level nears the band edges: consequently, large 
amounts of holes and electrons are generated on the op- 
posite sides of the sample. These intrinsic carriers screen 
partially the polarization charges, preventing the gap 
from closing. The total potential drop is thus pinned 
at the gap value for all thicknesses d > d c - that is, the 
effective gap decreases down to zero, but not below. For 
d > d c , the field will decrease as 

|E| = £gapM 

For this picture to hold, the spatial extension of the 
screening charge at the sample surface must be compa- 
rable_with that of the polarization charge (a few A at 
mosto) and much smaller that the sample size. This will 
cause the field inside the sample to remain uniform, since 
the net effect of screening will be to change the effective 



polarization charge. Indeed, this assumption turns out to 
be verified in practice on direct inspection, as we discuss 
below. 

Clearly, the above mechanism will strongly influence 
QW's of thicknesses equal to, or larger than, the critical 
value d c . For the system we are considering here, with a 
built-in field of -2.26 MV/cm, the critical value is d c ~ 
165 A. To confirm our picture, we simulated QW's with 
the same composition and geometry considered in Sub- 
sec. B, and thicknesses below and above d c , to mimic the 
crossover from a "microscopic" to a "macroscopic" sam- 
ple. In this case, we need to describe very extended bulk 
regions on the left and right of the QW, in order to ac- 
count for the large screening length. Thus, we have made 
use of a classical Thomas-Fermi model, where the charge 
densities are calculated with Fermi-Dirac statistics of a 
classical system rather than by solving the Schrodinger 
equation in the TB basis. This allows to consider devices 
with a spatial extension of several hundreds of microns. 
Effective masses in this calculation were fitted to accu- 
rately reproduce the self-consistent TB results. 

The resulting self-consistent potential is shown in Fig- 
ure | for well thicknesses of 100, 200, 300, and 400 A. A 
first point to note is that the field remains uniform for all 
well thicknesses. The field value equals the polarization 
field for the smallest thickness (smaller than d c ); for the 
thicker wells, the field (while remaining uniform) indeed 
decreases progressively as oc l/d. 
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FIG. 8. Field self-screening in a thick layer near and above 
the critical thickness. 

Photoluminescence experiments are not expected to be 
able to reveal this effect (which should cause a saturation 
of the red shift as function of thickness) in very thick 
QW's, since the effective recombination rate rapidly be- 
comes vanishingly small. Experiments aiming to reveal 
this effect should be designed considering our result, that 
a very thick layer is effectively subjected to a uniform 
electrostatic field E gSjV /d. In an unstrained GaN QW, 
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for d > d c (the latter being typically of order 100-200 A 
or so depending on the polarization) the field is ~3.4 V/d, 
i.e. ~70 kV/cm for d =0.5 ^m. This is presumably suf- 
ficient to cause observable bulk-like effects such as shifts 
in response functions or field effects on impurities. 

A similar "self-screening" behavior has been revealed 
indirectly in devices comprising sufficiently thick layers. 
In Ref. [52| a 300 A thick Alo.15Gao.85N layer was grown 
on a very thick GaN substrate, and topped with a Schot- 
tky contact. The predicted field in the AlGaN layer is 
1.4 MV/cm, which would cause a potential drop of 4.2 
eV across the layer. The maximum reasonable poten- 
tial drop dictated by Schottky barriers, conduction off- 
set, and Fermi level is about 1 eV, so it must be the 
case that the polarization charge gets largely screened 
by electrons from the GaN layer, forming a high-density 
two-dimensional electron gas (2DEG) at the heteroin- 
terface; this, by the way, causes an enhanced conduc- 



epth profiling indeed 
An equivalent, more 



tivity in the active channel. CV r 
reveals a 2DEG at the interfaced 
formal description is that the field would force the metal- 
determined Fermi level to some 3 eV above the conduc- 
tion band of GaN, thus attracting towards the interface 
an enormous carrier density, which screen out (part of) 
the field. Note in passing that in Ref. [32] only piezo- 
electric polarization was considered, which leads to an 
underestimation of the 2DEG density, since the piezo- 
electric contribution is actually about one third of the 
total interface charge. Siaailar considerations apply to 
other similar experiments. E3 The recent device simula- 
tion of Ref. |^ has corrected this point, including in part 
the spontaneous-polarization interface charges. 



D. Suppressing QCSE by doping 

We have seen in the previous Sections that polarization 
fields can be screened to a reasonable extent by genera- 
tion of free charge of both kinds in the QW upon e.g. op- 
tical excitation. Qualitative problems with this screening 
mechanism are that (a) it is transient, since it disappears 
when photoexcitation or current injection are removed, 
and that (b) in purely electronic {i.e. non-optoelectronic) 
devices, it is unlikely that the high densities needed can 
be reached in normal operating conditions. Besides, the 
current is not constant in time, so that the well shape 
also changes in time. 

It is natural to presume that the same effects can be 
achieved in a permanent fashion using extrinsic carriers 
from dopants. The idea is to provide the well with carri- 
ers which would screen the polarization charge, excepts 
that now the electrons are released into the QW from 
the doped barriers, and not injected or photogenerated. 
Of course, this effect is not transient as the others dis- 
cussed previously. The problem is, how high must the 
doping density be to achieve the same level of screening 
as in a high optical excitation regime. We simulated a 50 



A thick Alo.2Gao.8N/GaN single QW, where the barriers 
have been doped n- type in the range-ironi 10 17 to 10 20 
cm~ 3 and the donor ionization energyea have been set to 
10 meV. In this simulation, we used again the sclfconsis- 
tent TB approach. The resulting conduction band profile 
is displayed in Figure ^| for the various doping densities. 
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FIG. 9. Conduction band edge of a remotely-doped 
GaN/ AlGaN quantum well. 

The polarization field raises the conduction band on 
the left side over the Fermi energy, and in order for the 
barrier conduction band to reach the Fermi level on the 
far left, the electrons are transferred from the left-side 
barrier into the QW, leaving behind a large depletion 
layer. As a consequence of the electron flow into the QW, 
the polarization field start to get significantly screened at 
doping densities above ~ 10 19 cm~ 3 . The existence of a 
depletion layer causes a large band bending in the left- 
side barrier, while the bending is absent in the left half 
of the well. This is quite different to the case of the pho- 
toexcited well, where the bending on the left side of the 
well was due to hole accumulation near the interface (see 
Figure 6). This explains why the field remains nearly uni- 
form in the left half of the well for all of the simulations 
performed. On the right side of the well, only a small 
bending due to the electron accumulation is present. In- 
deed electron localization is quite weak in these systems, 
since the confinement potential is small as compared to 
the field-induced drop, and electrons tend to spill over 
to the right-side barrier. This is likely to be case in all 
nitride systems in this composition range. 

From these results, we conclude that doping can indeed 
be used to screen polarization fields. While it is not ob- 
vious that the needed doping level can always be reached 
in practice, it is likely that a combination of doping and 
current injection or photoexcitation will generally succed 
in quenching polarization fields in the range of MV/cm, 
thus allowing for recovery of quasi flat band conditions. 
Fields in InGaN/GaN systems will be generally smaller 







than those in AlGaN/GaN systems for typical compo- 
sitions in use today, and will therefore be more easily 
amenable to treatment by the above technique. This pro- 
cedure has in fact been adopted in experiment by Naka- 
mura's group,H2l which reported that a doping level of 
10 19 cm~ 3 is sufficient to quench the QCSE to a large 
extent. Indeed, in their Ino.15Gao.85N/GaN MQWs the 
unscreened field is ~1.2 MV/cm, i.e. approximately a 
half of the one we considered here. This field will be 
more easily screened by remote doping, in qualitative 
agreement with our findings. 



VI. SUMMARY AND ACKNOWLEDGEMENTS 

In conclusion, we have discussed how macroscopic (and 
in particular, spontaneous) polarization plays an impor- 
tant role in nitride-based MQWs by producing large 
built-in electric fields. Contrary to zincblende semicon- 
ductors, in III-V nitrides-based devices the spontaneous 
polarization is an unavoidable source of large electric 
fields even in lattice-matched (unstrained) systems. The 
existence of these fields may also be used as additional 
degree of freedom in device design: for instance, for an 
appropriate choice of alloy composition, spontaneous and 
piezoelectric fields may be caused to cancel out, thus free- 
ing the structures from built-in fields. We have also dis- 
cussed the different regimes of free carrier screening, ef- 
fected by doping or optical excitation, showing that fields 
can be screened only in the presence of high free carrier 
densities, which leads to unusually high lasing thresholds 
for undoped QW's. Of course, our results about the ef- 
fects on the electronic structure apply qualitatively to 
any kind of polarization field, thus in particular also to 
piezo-generated ones. 
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